Parameter-free expression for superconducting T c in cuprates 
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A parameter-free expression for the superconducting critical temperature of layered cuprates is 
derived which allows us to express T c in terms of experimentally measured parameters. It yields T c 
values observed in about 30 lanthanum, yttrium and mercury-based samples for different levels of 
doping. This remarkable agreement with the experiment as well as the unusual critical behaviour 
and the normal-state gap indicate that many cuprates are close to the Bose-Einstein condensation 
regime. 



1 



An ultimate goal of the theory of superconductivity is 
to provide an expression for T c as a function of some well- 
defined parameters characterizing the material |Q. The 
BCS theory provides a 'material' aspect in an estimate 
of T c since the characteristic phonon frequency and the 
electron-phonon coupling constant can be measured || 
while the Coulomb pseudopotcntial is determined in the 
first principle (LDA) band-structure calculations. There 
is no general restrictions on the BCS value of T c if the di- 
electric function formalism is properly applied (^^]. On 
the other hand, there is a growing evidence for a non- 
Fermi-liquid normal state and for a non-BCS supercon- 
ducting state of doped cuprates (see, for example ||). 
In particular, the correlation of T c with the carrier den- 
sity and with their effective mass HQ, the carrier spe- 
cific heat near the transition and its magnetic field de- 
pendence ||, the divergent upper critical field iPHllf are 
reminicent of a charged Bose liquid |||,|l2],[l3) ■ 

In this letter we derive a parameter-free expression for 
T c of a charged Bose-liquid on a quasi-2D lattice. It 
yields T c values observed in many cuprates for any level 
of doping. The main point of our letter is largely in- 
dependent of the microscopical interpretation of charged 
bosons, which might be lattice and/or spin bipolarons 
||, or any other preformed pairs. 

In the framework of this rather general model T c is cal- 
culated from the density sum rule as the Bose-Einstein 
condensation temperature of 2e charged bosons on a lat- 
tice. Just before the discovery [Q we estimated T c as 
high as ~ 100K by applying our expression for the bipo- 
laron effective mass That was tested by Uemura 

[plpl with the conclusion that cuprates are neither BCS 
nor BEC superfluids but they are in a crossover region 
from one to another. The experimental T c has been 
found about 3 or more times below the BEC tempera- 
ture. 

We now calculate the critical temperature of a charged 
Bose-liquid taking carefully into account the microscopic 
band structure of bosons in layered cuprates. It appears 
that due to a multi-band structure bosons have internal 
quantum number that might give a two-fold degeneracy 
as derived by one of us for a particular case of a 'peroxy' 
bipolaron Jnj. We arrive at the parameter- free expres- 
sion for T c , which involves the in-plane, X a b and out- 
of-plane, A c magnetic field penetration depths and the 
normal state Hall coefficient Rh just above the transi- 
tion. It describes the experimental data for a few dozen 
different samples clearly indicating that many cuprates 
are in the BEC regime. 

Bound pairs in cuprates are intersite pairs j^] because 
of a strong on-site repulsion. In the case of the electron- 
phonon coupling this is confirmed by the numerical sim- 
ulations of ionic perovskite lattices (f^| , where two types 
of pairs were found, i.e a 'peroxy' bipolaron and an 'in- 
plane' one. The energy spectrum of the 'peroxy bipo- 
laron ( a bound state of the in-plane and apical holes) is 



a doublet due to two (x and y) oxygen orbitals elongated 
along the Cu02 planes ]l7j . The energy band minima 
are found at the Brillouin zone boundary, (±7r, 0) and 
(0, ±7r) owing to the opposite sign of the ppa and ppir 
oxygen hopping integrals. Near these points an effective 
mass approximation is applied with the following result 
for the x and y bipolaron bands: 
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where d is the interplane distance and t c /2 is the inter- 
plane hopping integral. The wave vectors corresponding 
to the energy minima belong to the stars with two prongs. 
Their group has only ID representations. It means that 
the spectrum is degenerate with respect to the number 
of prongs of the star. The spectrum Eq.(l) belongs to 
the star with two prongs, and, hence it is a two-fold de- 
generate. 

While the X-ray absorption experiments confirm an 
important role of apical holes in some cuprates (T^| there 
are multilayer compounds with inner layers without apex 
oxygens. The in-plane p-band hybridised with copper 
might have a higher energy than other bands and the 
first to be doped. In those cases the in-plane boson ( 
a bound state of two oxygen holes in the CuOi plane) 
is the ground state with the energy minimum, which 
might be found at the T point of the Brillouin zone. One 
can apply the method of invariants to derive its spec- 
trum. The space group of Lai- x Sr x CuO± in the tetrag- 
onal phase is D^ h . The wave functions at the T point 
transform as ID A\ gtU , A% g , u , Bi g ^ u , B2 g , u or 2D E g and 
E u representations of the point group D^h- For the case 
of Eg (basis functions xz and yz) and E u (basis func- 
tions x and y) the spectrum is expressed as the eigenval- 
ues of the 2x2 matrix. This matrix can be written as a 
linear combination of Pauli matrices, r,-. Taking into ac- 
count that To transforms as A\ g (basis functions k 2 + kz, 
k 2 z ), t x transforms as B\ g (basis functions k 2 — ky), r y 
transforms as Big (basis functions k x k y ), t z transforms 
as A2 g (basis functions s z ) we write the matrix of the 
Hamiltonian as H = [A{k 2 x + ky) + Bk 2 ] r + C{k 2 x - 
ky)T x +Dk x k y T y . The eigenvalues of this Hamiltonian are 

£k 2 = A{k 2 x + k 2 y ) + Bk 2 z ) ± ^C 2 {k 2 x ~k 2 y ) 2 +D 2 k 2 x k 2 y . 
Here A, B, C, and D are fenomenological constants that 
parametrize the effective mass tensor. As a result the 
in-plane spectrum is degenerate in T point as well if it 
belongs to a 2D representation. 

It should be pointed out that a low temperature phase 
has lower symmetry (space group D^f.)- It means that 
the degeneracy of the spectrum might be removed below 
the structural phase transition. Nevertheless, if the level 
splitting is less or of the order of T c we can apply a theory 
with a two-fold degenerate spectrum. On the other hand, 
if the degeneracy is removed the theoretical T c would be 
higher by 2 2//3 (see below). 
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Substituting the spectrum, Eq.(l) into the density sum 
rule, 

J2 [eMK/kBTj-iy 1 =n B (2) 



one readily obtains T c as 
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where the coefficient / ~ 1 is shown in Fig.l as a func- 
tion of the anisotropy, t c /k B T c , and m c = ?i 2 /|f c |c? 2 .This 
expression is rather ambiguous so far because the ef- 
fective mass tensor as well as the boson density n B 
are unknown and doping dependent due to the screen- 
ing of the interaction by free carriers and their local- 
isation by disorder. Fortunately, one can express the 
band-structure parameters through the in-plane, X a b — 
[m x m y /8Trn B e 2 (m x + my)} 1 ' 2 and out-of-plane penetra- 
tion depth, A c — [m c / lQirnse 2 } 1 ^ 2 (we take c = 1). The 
boson density is expressed through the in-plane Hall con- 
stant (above the transition) as p7|: 
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which leads to 
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with T c measured in Kelvin, eR B in cm and A in cm. 
The coefficient / is about unity in a very wide range of the 
anisotropy t c /k B T c > 0.01, Fig.l. As a result we arrive 
at a parameter-free expression, which unambiguously can 
tell us how far cuprates are from the BEC regime: 



T c ~ T C (3D) = 1.64 
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It does not contain the mass tensor explicitly. Hence, 
any other dispersion law would lead to about the same 
result for T c if the spectrum is two-fold degenerate. 
As an example, taking for simplicity C = D = 
in the in-plane bound state dispersion, one arrives at 
T c ~ 3.31(n B /2) 2/3 /(m 2 mc) 1/3 , R H = l/2en B , and 
A a 6 = [m ab /16irn B e 2 } 1/2 , A c = [m c /167m B e 2 ] 1/2 , and 
with Eq.(6) as well. Here m a b — fr 2 /2A and m c = h 2 /2B. 

We compare the theoretical expressions, Eq.(5,6) with 
the experimental T c of more than 30 different cuprates, 
for which both A a h and A c are measured along with Rh 
(see Table 1 and Fig. 2). The Hall coefficient has a strong 
temperature dependence above T c . Therefore, we use the 
experimental Hall 'constant' R B = Rh(T c + 0) just above 
the transition as prescribed by Eqs.(5,6). In a few cases 
(mercury compounds), where Hh(T c + 0) is unknown, 
we have taken the inverse chemical density of carriers 



(divided by e) as R#. Almost for all samples the theo- 
retical T c coincides with the experimental one within the 
experimental error bar for the penetration depth (about 
±10%). There are a few Zn doped YBCO samples, Fig.2, 
with the experimental critical temperature higher than 
the theoretical one. We believe that the degeneracy of 
the boson spectrum is removed by the random potential 
of Zn, so for these samples the theoretical T c is actually 
higher than Eq.(6) suggests. Multiplying the theoretical 
T c in Table 1 by 2 2 / 3 for three Zn-doped samples we ob- 
tain T c =73 K, 57 K, 41 K in good agreement with the 
experimental values T c =68 K, 55 K and 46 K, respec- 
tively. 

One can argue that cuprates belong to a 2D 'XY' uni- 
versality class with the Kosterlitz-Thouless critical tem- 
perature Tkt of preformed bosons (2(]]2l[] or Cooper 
pairs |^2| due to a large anisotropy. Would it be the case 
one could hardly discriminate Cooper pairs with respect 
to local pairs in cuprates by the use of their T c values. 
The Kosterlitz-Thouless temperature is expressed trough 
the in-plane penetration depth alone as p2| 
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It appears significantly higher than the experimental 
value in most cases (see Table 1). What is more cru- 
cial, however, is the fact that cuprates have the specific 
heat H of a 3D charged Bose gas Q . 

The boson-boson interaction might be rather strong 
leading to self-energy effects and to some renormalisa- 
tion of the effective mass tensor. It is important, that 
Eq.(6) does not contain the mass and, hence, is not af- 
fected by the interaction. Nevertheless, it is interest- 
ing to evaluate the effective mass tensor in terms of the 
penetration depth and the Hall constant. The in-planc 
and out-of-plane boson masses are presented in Table 2 
for a few samples of La2- a; Sr 2 ;Cu4 and YBa2Cu3C , 7_ ;r . 
The in-plane boson dispersion with C = D = has 
been applied. Theoretical estimates of the in-plane mass 
(~ 10m e i.e.about 5m e per hole (0||||Jl|||) fit well 
our empirical values, Table 2. There is an intersting op- 
posite tendency in the doping dependence of the effective 
mass of La2- 2; Sr 2 ;Cu04 and YBa2Cu307_a;. While the 
mass increases with doping in L^-zSr^CuO^, it slightly 
decreases in YBa 2 Cu307_ ;E . We believe that it is a re- 
sult of an interplay between an interaction responsible for 
the mass enhancement and disorder. We notice, however, 
that the absolute value of the effective mass in terms of 
the free electron mass, Table 2, does not describe the 
actual band mass renormalisation because the unrenor- 
malised (bare) band mass remains unknown. 

Many thermodynamic, magnetic and kinetic proper- 
ties of cuprates were understood with charged bosons on 
a lattice j5 26|. We admit, however, that one experi- 
mental fact might be sufficient to destroy any theory. 
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In particular, the single-particle spectral function seen 
by ARPES was interpreted by several authors as 
a Fermi liquid feature of the normal state incompatible 
with bipolarons. Most (but not all) of these measure- 
ments produced a large Fermi surface. This should evolve 
with doping as (1 — x) in a clear contradiction with low 
frequency kinetics and thermodynamics, which show an 
evolution proportional to x (x is the number of holes in- 
troduced by doping). Recently it has been established, 
however, that there is a normal state gap in ARPES and 
SIN tunnelling, existing well above T c irrespective to the 
doping level |2^-|2S|] . 'The 'Fermi surface' showed missing 
segments just near the points [ p8| where we expect the 
Bose-Einstein condensation |k| . A plausible explanation 
is that there are two liquids in cuprates, the normal Fermi 
liquid and the charged Bose- liquid (this mixture has been 
theoretically discussed a long time ago |3l]]). However, 
it is difficult to see how this scenario could explain the 
doping dependence of dc and ac conductivity as well as of 
the magnetic susceptibility and carrier specific heat. On 
the other hand, the single-particle spectral function of a 
pure 'bosonic' system has been recently derived by one 
of us (|^] . It describes the spectral features of tunnelling 
and photoemission in cuprates. Any single-particle spec- 
tral weight at the chemical potential appears in our model 
due to states localised by disorder inside the normal state 
gap |30|. The model is thus compatible with the doping 
evolution of thermodynamic and kinetic properties. 

In conclusion, we have shown that the experimental 
critical temperate of superconducting cuprates is not very 
different from the Bose-Einstein condensation tempera- 
ture of two-fold degenerate charged bosons on a lattice. 
Our empirical expression for T c describes the experi- 
mental data remarkably well with no parameters to fit. 
This possibility originates in two different energy scales 
in cuprates: a strong attractive interaction and a small 
bandwidth. Their difference allows us to 'integrate out' 
the interaction and express T c via static response func- 
tions. 
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Figure Captures 

Fig.l. Correction coefficient to the 3D Bose-Einstein 
condensation temperature as a function of the anisotropy. 

Fig. 2. Theoretical critical temperature compared with 
the experiment ( the theory is exact for samples on 
the straight line) for LaSrCuO compounds (squares) 
for Zn substituted YBa^Cui-^Zn^O^ (circles) for 
YBa2Cu 3 07_(5 (triangles) and for HgBa 2 Cu04 +< 5 (dia- 
monds). Experimental data for the London penetration 
depth are taken from T. Xiang et al, Int. J. Mod. Phys. 
B 12, 1007 (1998) and B. Janossy et al, Physica C, 181, 
51 (1991) in YBa 2 Cu 3 07-,5 and YBa 2 (Cui_ x Zn x ) 3 07; 
from V.G. Grebennik et al, 'Hypcrfine Interactions', 61, 
1093 (1990) and C. Panagopoulos (private communica- 
tion) in underdoped and overdoped La2- x Sr ;c Cu04, re- 
spectively, and from J. Hofer et al, Physica C, 297, 
103 (1998) in HgBa 2 Cu0 4+< 5. The Hall coefficient above 
T c is taken from A. Carrington et al, Phys. Rev. B 
48, 13051 (1993) and J. R. Cooper (private communi- 
cation) (YBa 2 Cu 3 07_5 and YBa 2 (Cui_ x Zn x ) 3 07) and 
from H.Y. Hwang et al, Phys. Rev. Lett. 72, 2636 
(1994) (La 2 _ x Sr x Cu0 4 ). 

TABLE I. Experimental data on T C (K), ab and c 
penetration depth(nm), Hall coefficient (10 -3 (cm 3 /C)), 
and calculated values of T c applying Eqs.(6,5,7) respec- 
tively (K)for La 2 - x Sr x CuOi (La), YBaCuO(x%Zn) {Zn), 
YBaiCuzO-r-x (Y) and HgBaiCuO^+x (Hg) compounds 



Compound 


T c 


A a f) 


A c 


Rh, 


T c 


T c 


Tkt 


Lo(0.2) 


36.2 


200 


2540 


0.8 


38 


41 


93 


La(0.22) 


27.5 


198 


2620 


0.62 


35 


36 


95 


La(0.24) 


20.0 


205 


2590 


0.55 


32 


32 


88 


Lo(0.15) 


37.0 


240 
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1.7 


33 


39 


65 


Lo(0.1) 


30.0 


320 


4160 


4.0 


25 


31 


36 


La(0.25) 


24.0 


280 
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0.52 


17 


19 


47 


Zn(0) 


92.5 
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111 


114 


172 


Zn(2) 


68.2 
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46 


50 


Zn(3) 


55.0 
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36 


38 


Zn(5) 


46.4 


370 
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1.2 


26 


26 


30 


Y(0.3) 


66.0 


210 


4530 


1.75 


31 


51 


77 


L(0.43) 


56.0 


290 


7170 


1.45 


14 


28 


40 


F(0.08) 


91.5 


186 


1240 


1.7 


87 


88 


98 


F(0.12) 


87.9 


186 
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75 


82 


97 
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83.7 
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115 
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78.2 


161 
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8.2 


43 


92 
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78.5 


200 


12600 


8.2 


28 


69 
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6.85 


56 
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229 
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145 
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79 


120 
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95.3 


165 
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61 
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197 
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4220 
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66 


105 


216 


ff S (0.101) 


93.4 


156 


3980 


4.48 


70 


107 


220 



#p(0.101) 92.5 139 3480 4.4 88 127 277 
#g(0.105) 90.9 156 3920 4.3 69 106 220 
##(0.108) 89.1 177 3980 4.2 58 90 171 



TABLE II. Mass enhancement with respect to the free elec- 
tron mass 

Compound m ab m c 

La(0.2) 22~J 3558 

La(0.15) 15.0 2698 

La(0.1) 11.3 1909 

Y(0.0) 7.2 584 

L(0.12) 8.3 600 

y(0.3) 106 1994 
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